Comment on "Mean First Passage Time for Anomalous Diffusion" 
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We correct a previously erroneous calculation [Phys. Rev. E 62, 6065 (2000)] of the mean first 
passage time of a subdiffusive process to reach either end of a finite interval in one dimension. The 
mean first passage time is in fact infinite. 
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Anomalous diffusion is commonly characterized by the 
behavior of the mean squared displacement as a function 
of time at long timesP,0,l3, 



(x 2 ) 



2K a 
T(l + a) 



(1) 



where K a is a generalized diffusion constant. Ordinary 
diffusion corresponds to a = 1 and, in the more usual 
notation, K\ = D. Superdiffusion is associated with mo- 
tion that is faster than diffusive, a > 1, while subdif- 
fusion occurs when a < 1. In a recent paper £j, the 
mean first passage times to the ends of an interval for 
a superdiffusive and a subdiffusive random walker on a 
line were calculated, and these results have even more 
recently been applied to the problem of anomalous heat 
conduction in such a line in the presence of a temperature 
gradient However, there is an error in the calculation 
for the subdiffusive problem from which one concludes 
that the mean first passage time reported in 4] for this 
case are incorrect. In fact, the mean first passage time 
for the subdiffusive problem is infinite. 

To support this observation it is compelling to note 
that a continuous time subdiffusive nearest neighbor ran- 
dom walk (CTRW) |3( with a waiting time distribution 
which has a long tail, i.e., a walk in which the probabil- 
ity density that a particle takes the next step at a time 
t — > oo after the previous step is ip(t) ~ C a /t 1+a leads 
exactly to Eq. JIJ (C Q is a constant). The mean time for 
the particle to make a single jump is 



where 



For large t, 



T\ = lim Ti(t), 



T 1 (t)= / r^(r)dr. 
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so that one obtains the well-known result T\ = oo, i.e., 
the mean time to go from any one location to another 



even in a single jump is infinite. However, this argument 
might generate issues about the waiting time for the first 
step of the process, since one of the differences between a 
CTRW and the fractional diffusion equation lies precisely 
in the assumptions associated with this first step. In a 
CTRW there is a singular contribution to the probability 
density that the particle is still at the origin xo at time 
t. 



P(x,t) 



t a S(x ) + other terms, 



(5) 



which does not appear in the solution of the fractional 
diffusion equation 

To sidestep this problem and show that the divergence 
of the mean first passage time does not arise only from 
this term, we also obtain the divergent result starting 
with the fractional diffusion equation that was the start- 
ing point in |4|. Although his general formulation was 
for an arbitrary starting site in the interval (0, L) and in 
the presence of an external constant force, the explicit 
final result was presented for a particular initial location, 
x = L/2, and with no external force. This explicit result 
is also the one used in 0. We thus restrict our presen- 
tation to this specific case. 

The mean first passage time from x = L/2 to either 
x = or x — L is given by 

/>L />oo />oo 

T= dx dtP(x,t)= / dtS(t), (6) 
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where 



S(t) = / dxP(x,t) 
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and P[xJ;)Js the solution of the fractional diffusion equa- 
tion 



K aa D]- a ^P{x,t) 



(8) 



with absorbing boundary conditions P(0,t) = P(L,t) = 
and initial condition P(x,t = 0) = S(x — L/2). Here 
Q n - Q i s the Riemann-Liouville operator 
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and K a is the generalized diffusion coefficient in Eq. . 
The quantity S(t) is called the survival probability be- 
cause, as one sees from Eq. |0, S(t) is just the prob- 
ability that the particle has not been absorbed by the 
boundaries at x — and x = L during the time interval 

[CM]. 

The solution of Eq. JHJ with the given boundary and 
initial conditions can be found by the method of separa- 
tion of variables 0, 0] : 
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Here E a (—z) is the Mittag-Leffler function (for a = 1 it 
reduces to the exponential exp(— z) and thus yields the 
usual solution for the diffusive problem). It then follows 
that 
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The mean first passage time to s = or L is then T 
limt^oo T(t) where 



T(t)= f AtS{t). 
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To address the convergence of T(t) for t — > oo we need 
to analyze the long-time behavior of S(t). Note that S(t) 
is well behaved for finite times (the survival probability 
goes to 1 for t — > 0), so that the divergence of T is due to 
the behavior at long times. For large z the Mittag-Leffler 
function behaves as 
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and consequently, for large t, 
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where Z(m) = T,n=o(- 1 ) n / ( 2n + l) 2m+1 . For f -> oo, 
and using the fact that 2(1) = 7r 3 /32, we then have 



S(t) 
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8r(l - a) X a t Q ' 
It then follows that for large t we have 



T(t) = f 
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i.e., T(t) oc exactly as in Eq. Q. We thus conclude 

that T(t) — > oo when t — * oo for any a < 1. 

We have thus shown that the mean first passage time 
for a subdiffusive process described by the fractional dif- 
fusion equation (or, for that matter, by a continuous 
time random walk) to reach the boundaries of a one- 
dimensional interval is infinite. We note that our Eq. i|ll|) 
appears as Eq. (40) in [TJ, but the connection between 
the survival probability and the mean first passage time 
is never made in that work so a user of the result in £| 
would not necessarily discover the connection. An ex- 
pression for the first passage time density involving the 
Mittag-Leffler function appears as Eq. (3.87) in [lfj, but, 
again, they do not go on to calculate the mean first pas- 
sage time, nor do they do the necessary asymptotic anal- 
ysis of the Mittag-Leffler function that would allow them 
to do so. That these results would not necessarily lead 
a reader to conclude that the mean first passage time is 
infinite is reinforced by the fact that both of these refer- 
ences appear in Q. 

Finally, we note that the validity of the result in |4| for 
the mean first passage time in the superdiffusive regime 
has also been questioned recentl y becau se it violates a 
theorem due to Sparre Andersen 
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